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Abstract 

Based on the Dirac equation, the behavior of relativistic electrons which tunnel a potential barrier 
of height Vo for incoming energies between Vo and Vo + m, is studied by using the wave packet 
formalism. The choice of this incoming energy zone guarantees that only electrons participate to 
the tunneling process. In the opaque limit, as shown in a previous analysis, the transit velocity is 
proportional to the barrier width and for relativistic potential tends to a constant value greater 
than c. In this paper, a new numerical study shows a very surprising result: superluminal 
transmissions are also evident for thin barriers. 

I. INTRODUCTION 

In this paper, by using the Dirac equation, we examine the relationship between the peak of the 
incoming wave packet (of energy E and mass m) which encounters a potential barrier of height Vo 
and width L at z ~ and the peak of transmitted wave packet which moves in the free potential 
region after the barrier. The terminology "Dirac tunneling" refers to the energy zone of the incoming 
particles, Vq < E < Vo + m [T]-[3]. This energy zone is characterized by evanescent solutions in the 
region < z < L, i.e. exp[±pz] with p = yni^ — {E — VoY . The terminology Dirac tunneling is 
useful to distinguish these evanescent solutions form the evanescent solutions which appear in the 
Klein tunneling zone, Vo — m < E < Vo [3]. The study of the Klein tunneling zone exceeds the scope 
of this paper and it will be appropriately discussed in a forthcoming article. It is important to observe 
that the Klein tunneling zone has not to be confused with the standard Klein zone, E < Vo — m, 
which is characterized by oscillating antiparticle solutions, i.e. exp[ztiqz] with q = \J {E — VoY — m^ 
and the appearance of the Klein paradox [4HZ]- For the convenience of the reader, we draw the energy 
Dirac/Klein zones in the following picture 
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The incoming electrons will be described by a wave packet with a gaussian momentum distribution, 
g{p) centered at po, 



dip) 



cxp - {p - Po) d /4 for y^Vf^ — ml 



otherwise , 



< p < p„ 



v/K(K+2TOe) 



where the minimum and maximum values of the momentum distribution are chosen to guarantee 
Dirac tunneling. The barrier filter effect modifies the incoming momentum distribution [SUllj. The 
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momentum mean value Po after transmission tends to an higher value, say pj^. In the opaque limit 
(jrieL 3> 1) pjj ~ p„„^. For thin barriers pj^ ss po- The incoming electrons move in the free region 
before the barrier with a subluminal velocity given by Va ~ po/Eg. The outgoing electrons move in 
the free region after the barrier with a greater, but obviously still subluminal, velocity vj^ = pj^/Ej^. 
This implies that the maximum of the spatial distribution of electrons, which tunnels a barrier of 
width L, will reach the position D > L at time 

L D- L 

''■^'-.^^^- « 

A natural question is to ask if and in which cases the average velocity v^^j — ^ I ^L ^^ superluminal. 
Our paper was intended as an attempt to give a satisfactory answer to this intriguing question not 
only in the opaque limit, where the transmission probability goes to zero, but also for thinner barriers, 
where greater transmission probabilities are found. Finally, the choice of the Dirac tunneling energy 
zone seems to be the more appropriate in discussing superluminal velocities because it avoids the 
possibility that the electrons which appear in the free potential region after the barrier are products 
of pairs production (a typical phenomenon of the Klein zone). The lower limit (% < rn) is not covered 
in this work. 

II. THE DIRAC EQUATION IN PRESENCE OF A POTENTIAL BARRIER 

Since we shall need the solutions for a barrier potential, we rewrite the free Dirac equation 

i 7''a^*(r, t) = me *(r, t) , (2) 

by including, via minimal coupling 9^ ^ 9^; — i e A^ p^l3| . an electrostatic potential A^ = \A{z) , 0], 
where A{z) — A^ for < z < L and otherwise, 

[ 1 7^a^ + e 7o A(z) ] *(r , i) = me *(r , t) . (3) 

Considering a one-dimensional motion, px = Py = and pz ~ p, and stationary state solutions, 
ip{z)exp[—iEt)], we obtain 

-t„[E- Viz) ] i,{z) + 1 73 ^j'{z) = me ^(z) , (4) 

where V{z) = ~eA{z) and the prime indicates the derivative with respect to z. In any constant 
potential region, for a given E, only two solutions exist be they oscillatory (|ii^— Vo| > m) or evanescent 
(|i? — VIjI < m). Using the Pauli-Dirac set of gamma matrices and observing that for one-dimensional 
phenomena spin flip does not occur, we find, for Dirac tunneling, the following spinorial solutions [3] 

Region I : u{p, E) e+'P"" + Ru{-p, E) c-'p^ , 

Region II : Au{ip, E ~Vo) e'P'' + B u{-ip, E -Vo) e+P"" , (5) 

Region III : T u{p , S) e+'P^ , 

where 

u(p,S) = [l,0,p/(S + 7«e),0]* . 

Because the Dirac equation is a first order equation in the spatial derivatives, for a step- wise conti- 
nuous potential, only the continuity of the wave function has to be required. Imposing the continuity 
conditions, V'i(O) == V'i/(0) and ipi,{L) ~ iJji,i{L), we obtain the following transmission coefficient 



Tip,L) = e-'P'^/ 



., ,. . P' - VqE . ,, ^. 

cosh( pL) — I sinh( pL) 

pp 



(6) 



The transmitted wave packet is then obtained by integrating over all the possible stationary states 
modulated by the weighting function .g(p), 

1',(z,i) = r^dpg{p)T{p,L)u{p,E)e'iP'-''*^ 



Pr, 



dp[gAp,L),0,Mp,L),0f ,:0 ^^lpiz^L)^Et] ^ (7) 



where 



tan 6* 



p' - VqE 
PR 



tanh(pi) , g^{p,L)^g{p)\T{p,L)\ and fr{p,L) = 



P 



E + m. 



9t{p,L) . 



III. THE FILTER EFFECT FOR OPAQUE BARRIERS 

This section contains a brief summary of the discussion recently appeared in hterature [14] on the filter 
effect for opaque barriers in the Dirac equation. The action of the filter effect [5] on the momentum 
distribution of the transmitted wave, grip, L) and /r(p, i), is shown in Fig. 1. For increasing values of 
the barrier width L, the transmitted momentum averaged over the distributions Qt^J), L) and friPi L), 



PL •= (P>^ 



dpp 



glip, L) + U{p, L) 



dp 



glip, L) + f^{p, L) 



(8) 



tends to p^a^- This splitting of the average momentum is due to the higher momenta selection caused 
by the exponential functional dependence, exp[— pL], of the transmission coefficient T(p, L). The plots 
in Fig. 1 refer to a potential barrier of height Vf, — rrig, which implies p„i„ — and p„„^ = VS^^ej 
and an incoming momentum po = v3'Tie/2. The space localization for the incoming wave packet is 
determined by Tried = 10. The ratio between the down and up distribution, as a consequence of the 
higher momenta filter, increases from po/{Eo + rag) sa 37.3% to p^a^ / {Emax + "me) ~ 57.7%. 

Without going into the details of the mathematical apparatus of asymptotic integral expansions 
and, in particular, of the validity and possible generalization of the stationary phase method which 
find a complete description in a number of books }15 p i6 j and papers [8Ull| , wc briefly discuss the case 
of opaque barriers, i.e. rUeL ^ 1. In this limit, the cosh and sinh functions with argument pL can be 
approximated by exp[— pL]/2. This allows to rewrite the transmitted coefficient as 



\T{p,meL > 1)1 « 2ppe-P^ /meVo 



(9) 



Now, the computation of the integral given in Eq. ([7]), taking into account the previous approximation, 
can be done analytically. Starting from 



\'^T{L,t;meL:^ 1) 



gip^a.) 



2E„ 



rrieVo ' V E„ 

changing the variable of integration from p to p, 



dppp e-''^e^(«-^*) 



(10) 



pdp ; 



E„ 



and using the following scries expansions for 9 and E, 

« 6'„„, - a,p + 0[p'] and E = E, 



pdp, p{p,^,^)^me and p(p„„,) = 



02 p- 



0[P1 



where a^ ~ \JVo{Vo + 2 Trig) /meVf, and Oa = l/2me, wc find 



|*T(i,i;meL> 1) 



dp p^ e '' exp [ i ( 02^ p^ — fli p )] 



dp p^ e 



2 p- pL 



1 + i ( 02^ P^ ^ Oi p ) — ( 02^ P^ ^ 0,1 p) / 2 



(11) 



It is clear from this discussion that the problem of finding the maximum of the transmitted wave 
packet reduces, in the opaque barrier limit, to find the maximum of the function 



■PrW = I s(2) - [(02^)' s(6) + a\s{A)] / 2 + 0102^ .s(5) + i [a^t s(4) - a^s{2,)] 



(12) 



where 



s{n) 



dpp"e 



" „- pi 



n\/ L 



The condition dt\\'^T{L,t;7JieL ^ 1)|| = is thus satisfied requiring 7'^(i) — 0. Observing that 

Vrit) « s'{2) + al [s^(3) - s(2)s(4)] +2aia,i[s(2)s(5) - s(3)s(4)] + (a,t)^[s^(4) - s(2)s(6)] , 
we find for the tunneling time, Tt„„, the following equation 

ai[s(2)s(5) - s(3)s(4)] + a,T,^„ [s^(4) - s(2)s(6)] = . (13) 

4-4!/L^ -36-41. /L^° 

From this equation, we obtain a closed formula for the tunneling velocity. 



^ _Z^ ^ 9a, ^ 9 / Vo 



In order to understand the importance of Eq. (J14p . let us present two limit cases. For non relativistic 
potentials, 

l'*„„[K « ™e] « TT \hr^ « 1 . 

2 V 2me 
This means that the use the Schrodinger equation leads to non-relativistic tunneling velocities. A 
surprising result is seen in the chiral limit, 

9 
where superluminal tunneling velocities appear |14j . 

IV. NUMERICAL ANALYSIS AND TRANSIT TIMES 

Up to this point, we have limited ourselves to a strictly analytic approach to the problem. To test 
our analysis, we have performed a numerical study by using Mathematica 8 [Wolfram] whose results 
are shown in Fig. 2 and Fig. 3. The numerical analysis shows the phenomenon of multiple peaks. 
The spatial distribution of these peaks is plotted in Fig. 2 for a barrier of height Vo = rUe and 
incoming momentum distribution centered to po = v3r7T,e/2. In Table 1, we give the values of the 
transmitted probability density and the corresponding appearance time for the first peaks around 
the main maximum. For thinner barrier the secondary peaks are very small when compared to the 
main peak. The phenomenon of multiple peaks disappears for increasing values of the barrier width. 
In Fig. 3, we show the barrier width dependence of the tunneling time, Tt„„, referring to the main 
peak and the corresponding tunneling velocity, «(„„. The analytic value of the tunneling velocity, 
calculated by using Eq. p4|) for Vo = rrie, 9/2-\/3, is in good agreement with the numerical data. The 
important point to be noted here is that the numerical data clearly show superluminal tunneling 
velocities also for thinner barriers, see Fig. 3 in the initial range of the barrier width. The conclusion 
of our numerical study is unexpected and quite surprising. For thinner barrier superluminal tunneling 
velocities appear and the superluminal effects are greater than those observed in the opaque barrier 
limit. Due to this surprising result, we have investigated in details the behavior of an incoming wave 
packet of momentum Po = \/3me/ 2 and localization d = IO/ttIe which starting from the point z = at 
time t = reaches the point D = 40/me travelling through free space (Fig. 4a), tunneling a barrier of 
width L = lO/we (Fig. 4b), a barrier of width L = 20/me (Fig. 4c), and a barrier of width L = 30/TOe 
(Fig. 4d). The transmitted peak is found a.t D ~ iO/nie at the following times, {t^ j^}, 

{ rne tio,o , me ^40,10 , "le ^40,20 , me ^40, 30 } = { 61.35 , 46.03 , 29.95 , 15.70 } . 
These numerical times imply the following transit velocities, {v j. }, 

{ W40,0 , ^40,10 , ^40,20 , ^40,30 } == { 0.65 , 0.87 , 1.34 , 2.55 } . 

The condition for superluminality in D is 

vn, L = / ^, > 1 ^ D < — -^ i^ L . 



V. CONCLUSIONS 

The discussion of the time spent by a particle to tunnel through a barrier potential has been ap- 
proached from many different points of view and has been the issue of an intriguing debate in the last 
decades [TUJ[TT]. Since the energy is bounded from below, the time operator [T7] is not self-adjoint, 
there is no general agreement about a satisfactory definition of tunneling times and a universal intrinsic 
tunneling time valid for all experiments probably does not exist. 

The phase time is the time between the instant in which the peak of the incoming wave packet 
reaches the barrier and the instant in which the peak of the transmitted wave packet appears in the 
free region after the barrier. By using the stationary phase method, Hartman [5] observed that the 
tunneling thorugh opaque barriers is independent of the barrier length. This surprising result is due 
to the fact that the stationary phase method was applied without considering the filter effect. In the 
study presented in ref. |14| and briefly summarized in section III, the integral of the transmitted wave 
function for Dirac particles was calculated taking into account the filter effect. The new formula for 
the tunneling velocity in presence of opaque barriers. 



op.6... 9 / Vo 



V = — \ T c , 

2 V K + 2 m^c ' 

clearly shows that for non-relativistic potentials, Vq <^ rrieC^, the tunneling velocities are subluminal. 
There is not Hartman effect for the Schrodinger equation. Nevertheless, "superluminality" for the 
phase time is seen for relativistic potentials. For example, for Vg = m^c^ the tunneling velocity, given 
by 9 c / 2v3 , is superluminal and this is confirmed by numerical calculations, sec Fig. 2. The numerical 
data also present a still more surprising result. In Fig 3, it is clear that superluminal tunneling 
velocities appear not only in the opaque limit (very wide barriers) but also for thinner barriers. This 
stimulated a numerical analysis in which the transmission across free space, see Fig. 4. a, and through 
potential barriers of different widths, see Fig. 4.b-c, is compared. Superluminal transmissions appear. 
This apparent paradox clearly shows that the phase time is not the correct tunneling time definition. 
Why? The peak of the transmitted wave packet could be not related to the peak of the incoming 
one. When a wave packet with a given momentum distribution strikes a barrier, the transmitted wave 
packet will exhibit a distribution displaced to higher momenta. Consequently, the transmitted packet 
moves faster that the incident packet [T3] . Observing that the transmitted coefficient only depends on 
the barrier width and not on its position, see Eq. ^ , the transit time does not change if we displace 
the barrier from (0,L) to {a, a + L). This means, for example, that the time in which the peak of 
the transmitted wave packet reaches the point z = 40 h/m^ c in Fig. 4 is the same for barriers located 
at z = or at z > 0. Due to the fact that for z > the free region after the barrier is reduced, 
the outgoing peak cannot be related to the incoming one. Once again the filter effect explains the 
apparent paradox. Only higher momenta in the incoming packet contribute to tunneling. 

In view of the last comment, it could be interesting to review the results reported in this paper 
in terms of other definitions of tunneling times. The dwell time [18] is defined as the ratio of the 
number of particles within the barrier to the incident flux. The average of the time spent by the 
particle inside the barrier obviously does not distinguish whether, at the end of the process, the 
particles arc reflected or transmitted. Due to the fact that the effects shown in this paper imply a 
very small transmission, the average dwell time should be related, in this case, to the well-known 
delay time of quantum mechanics |15| (obviously using the Dirac formalism it should contain the 
appropriated relativistic corrections). A more complicated analysis is required if we introduce a time- 
modulated barrier |19| . V{t) = K) -|- VJ, coswi. In non-relativistic quantum mechanics, the crossover 
between small and large frequencies yields to the traversal time. The study of emitted and absorbed 
modulation quanta should be investigated in terms of relativistic transmitted particles. In this spirit, 
it could be also interesting to consider relativistic corrections to Larmor precession when x-polarizcd 
particles which move in the z-direction encounter in the barrier a small magnetic flcld pointing in 
the y-direction. The precession angle determines the time spent by the particle in the barrier [19]. 
Nevertheless, in recent papers, it was observed, for relativistic particles, the phenomenon of spin [20] 
and helicity |21] flip for planar motion. Thus, it could be very interesting to investigate, for Dirac 
planar motion, possible consequences in the spin and helicity flip related to transversal times. We 



conclude this discussion with our last suggestions on possible future studies in which the numerical 
analysis done in this paper could find possible applications. In recent papers |221I23| . the tunneling 
time problem for non-relativistic rectangular potential barriers was analyzed in terms of the tempus 
operator |24[|25j , —ihdj^. The tunneling time is expressed by an energy average of the phase time 
and there is no fundamental problem with the Hartman effect [22| . This conclusion is similar to our 
analysis for non-relativistic potentials. In [53], it was shown that when the energy of the incoming 
particle is near to a resonance pole of the tunneling barrier, an intrinsic tunneling time does exist, but 
when the energy is near the branch point there is no intrinsic tunneling time. It could be interesting 
to calculate the expectation value of the time operator, the resonance poles and the branch points 
for relativistic particles. Observing that the local value of the tempus operator gives a complex time 
for a particle to traverse a barrier, it should be also interesting to study the connection between 
the results obtained in the presence time formalism and the results obtained by introducing complex 
potential [^[?7] . In particular with respect to relativistic dccohcrence. 

We hope that the results presented in this paper could stimulate and motivate studies in tunneling 
phenomena by using the Dirac equation. In this spirit, this work can be seen as an initial step in view 
of new and more detailed discussions on the fascinating topic of relativistic tunneling. 
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Table 1: Transmission probability densities calculated at z — L and corresponding times in which the 
transmitted peaks appear in z = L. For thin barriers, the secondary peaks are very small with respect 
to the central peak. The phenomenon of multiple peaks disappears for increasing values of the barrier 
width. The numerical data of this Tabic was prepared by using an incoming momentum distribution 
centered at po = \/3mf,/2, a localization determined hy d = lO/nie, and a barrier potential of height 
Vo = me- 
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Figure 1: Momentum distributions for the up, g^iPjL), and down, fj.{p,L), component of the trans- 
mitted Dirac spinors. The incoming momentum distribution g{p) is centered at po = v3'Tie/2. Its 
space locahzation is determined hy d = lO/rrif.. By increasing the value of the barrier width L, the 
filter effect modifies the initial momentum distributions leading, for rrieL ^ 1, to smaller transmission 
probabilities, to greater momentum mean values, p-^ — >■ p^a.^, and to localizations ^ d. 
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Figure 2: Maximum and m.inimum. values of the transmitted probability density at z := L in correspon- 
dence to the time in which it appears in the free region after the barrier. The black circles indicates 
the central peak. The secondary peaks are identified by the black squares. Finally, the white squares 
represent the times in which the minima appear in z ~ L. The numerical data refer to incoming 
electrons, with a momentum distribution centered at Po = \/3me/2 and localization determined by 
d = lO/nie, which tunnel a potential barrier of height Vq = rUe and width L. 
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Figure 3: Tunneling times (black circles and continuous line) and tunneling velocity (withe circles and 
dotted line) for incoming electrons of momentum po = \/3me/2 and localization d = lO/irie which 
tunnel a potential barrier of height Vo ~ rrif, and width L. The numerical data confirm the analytic 
formula obtained in the opaque limit for the tunneling velocity, see Eq. (|14p . An unexpected and 
surprising result is seen for thin barriers (« 2 d) where the superluminal tunneling velocity reaches its 
maximum value (« 10). 
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Figure 4: Probability density for incoming electrons, with a momentum distribution centered at 
Po = \imf^/2 and localization d = 10/?7ie, which move in free space (a), tunnel a potential barrier 
of width L equal to lO/me (b), 20/ m^ (c), and 30/me (d). In the free space the electrons move 
with a mean velocity ^40,0 ~ 0.65. Barrier tunneling increases the velocity of propagation leading to 
superluminal transmission, ^40^20 ~ 1-34 and ^40^30 « 2.55. 



